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Abstract Magnetic fields appear to be ubiquitous in astrophysical environ- 
ments. Their existence in the intracluster medium is established through 
0^ ' observations of synchrotron emission and Faraday rotation. On the other 

hand, the nature of magnetic fields outside of clusters, where observations 
are scarce and controversial, remains largely unknown. In this chapter, we 
review recent developments in our understanding of the nature and origin of 
intergalactic magnetic fields, and in particular, intercluster fields. A plausi- 
ble scenario for the origin of galactic and intergalactic magnetic fields is for 
(~| , seed fields, created in the early universe, to be amplified by turbulent fiows 

C^H' induced during the formation of the large scale structure. We present several 

mechanisms for the generation of seed fields both before and after recom- 
. bination. We then discuss the evolution and role of magnetic fields during 

^ ' the formation of the first starts. We describe the turbulent amplification of 

, seed fields during the formation of large scale structure and the nature of the 

magnetic fields that arise. Finally, we discuss implications of intergalactic 
magnetic fields. 
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1 Introduction 



In the highly successful ylCDM cosmology, the large scale structure (LSS) 
forms through a process known as hierarchical clustering in which small-scale 
objects collapse first and merge to form systems of ever-increasing size. This 
scenario leads to a "cosmic web" of structure where galaxies reside mainly 
along filaments wh ile galaxy clusters arise at the intersections of two or more 
filaments (see, e.g. jBond et aLlflOO d). Furthermore, a picture of a multi-phase 
intergalactic medium (IGM) has emerged. A hot phase, often referred to as 
the intracluster medium (ICM) because it is found inside and around clusters 
and groups of galaxies, has a temperature T > 10^ K and is observable 
via X-ray emission. The warm-hot intergalactic m edium (WHIM) (10^ < 
T < 10^ K) is fou nd mainly in filaments of galaxies ()Cen and Ostrike"nll999l 
iKang et"aL 1 1200.^ . Gas in the IGM i s heated by cosmolo gical shocks which 
arise as objects in the hierarchy form. ( Rvu et aL ||20 3 , P fr ommer et al.l[2006L 



iKang et al.lboOTl ISkillman et al]|2008llVazza et al.ll2009il 



Magnetic fields are observed in galaxies of all types and in galaxy clus- 
ters. Moreover, there is some evidence that they permeate the filaments of 
the cosmic web. Various scenarios for the origin of galactic magnetic fields are 
discussed in IWidrow et ai] !i201Ct) of this volume. The basic idea is that weak 



fields, generated either by an exotic early universe mechanism or some astro- 
physical process, are amplified to fiG strength during galaxy formation and 
by dynamo action during the subsequent quiescent phase of galaxy evolution. 
Similarly, weak seed magnetic fields can be amplified into the intergalactic 
magnetic field (IGMF) by turbulent flow motions during the formation of the 
LSS through a process known as small-scale turbulence dynamo. In a tur- 
bulence dynamo, kinetic energy of the fluid in converted to m agnetic energy 
through stretching, twist ing, and folding of th e field (see, e.g.. ^Sub ramanian 



1999, Cho and Vislmia '3 l2000l iHaugen et al.l 2004a.b.c. Schckoch ihin et al 



2004I.IBrandenburg and Subramanian"2005l Cho et al.l2009llBrandenburg et al.l 



20101 of this volume). In addition, feedback from the black hole regions in 



active galactic nuclei (AGNs) can also contribute and amplify magnetic fields 
in the IGM. 

This chapter will explore the structure formation-magnetic field connec- 
tion. In particular, we will address two interrelated questions. First, can the 
structure formation itself generate new magnetic fields and amplify existing 
ones? Second, what role do magnetic fields play in structure formation and 
in other astrophysical processes? An outline of the chapter is as follows: In 
Section 2, we briefly review the observational evidence for magnetic fields 
in the clusters and filaments of the cosmic web. In Section 3, we describe 
various astrophysical mechanisms for the generation of seed magnetic fields, 
which would later be maintained and amplified within galaxies and clusters. 
In Section 4, we discuss the evolution and role of magnetic fields during the 
first star formation: their stabilization due to a non-zero ionization degree 
as well as their potential implication for star formation. We then describe 
how magnetic fields are amplified during structure formation in Section 5. 
We discuss the effects of magnetic fields on the propagation of cosmic rays in 
the IGM and also the Faraday rotation measure (RM) induced by magnetic 
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fields in Section 6. Section 7 presents something of a departure from tlie more 
phenomenological and astrophysical discussion as it provides a formal treat- 
ment of the evolution of density perturbations in the presence of magnetic 
fields. Finally, a brief summary follows in Section 8. 



2 Observational Evidence for Magnetic Fields in Clusters and 
Filaments 

The existence and strength of the IGMF has been deduced from observa- 
tions of s ynchrotron emission and confirmed directly from ob servations of RM 
(see, e-g.- ICarilh and Tavloill200l iGovoni and Ferettill2004l for r eview). Syn- 
chrotron emission from a galaxy cluster was first discovered by iLarge et al.l 
((1959) who surveyed the Coma cluster in the radio. Since then, it has been 
observed in numerous clusters , either as radi o halos or as radio relics (see, 
e.g., IGovoni and Ferettil 12004 [ Cassano et al.l[2008l . and references therein). 
The diffuse radiation from radio halos is mostly unpolarized, and is thought 
to emerge from turbulent magnetic fields in the ICM. The radiation from 
radio relics is highly polarized and is believed to be emitted from shocked 
regions in the ICM. The strength of the magnetic field in radio halos is es- 
timated to be ^ 1/iG, while stronger fields are found radio relics. These 
estimates either assume equipartition (also known as the minimum energy 
argument) or incorporate measurements of inverse- Compton emission. 

Observations of the IGMF based on Faraday rota tion have also been done , 
though mostly for magnetic fields in the ICM (see ICarilh and Tavloill2002L 
and references therein). An RM map of the Coma cluster, for example, re- 
veals a field with a strength ~ fiG and a coherence length of order ^ 10 kpc 
(|Kim et al.lll990l) . For Abeh clusters, the typical RM is ~ 100-200 rad m"^ 



which indicates an average field strength of ~ 5-10 fiG ("Clar ke et al.ll20dlL 
[Ciarke„2004) . RM maps of clusters can be used to determine the power spec- 
trum of turbulent magnetic fields in the ICM. For instance, a Kolmogorov- 
like spectrum with a bendin g at a few kp c scale is found in the cooled core 



region of the Hydra cluster (jVogt and E nfilin 2005) , and spectra consistent 



with the Kol mogorov spectrum w ere reported in the wider IC M for the Abell 
2382 cluster ([Guidetti et al.|[2008l ) and for the Coma cluster (jBonafede et al.l 
l201flh . 

The nature of the IGMF in filaments, on the contrary, remains largely 
unknown, because the studies of synchrotron emission and RM outside clus- 
ters are still scarce and controversial. Althoug h faint radio emi ssion has been 
observed in the outskirts of clusters (see, e.g.. lKim et aLllToSOl ). there are no 
confirmed observations of synchrotron emission from filaments. Such mea- 
surements present a challenge for current facilities. As well, the remov al of 
the galactic foreground is a non-trivial task (see, e.g., [Brown et al.ll2010h . At 
present there is an upper limit of ~ 0.1 /iG for the strength of the IGMF 
i n filaments based on the observ ed limit of the RMs of background quasars 
(|Rvu et al.l l l998llXu et al.ll200e 



Recentlv JNeronov and Vovkl (|2010t ) , lAleksic et all (|2010f ) reported a lower 



bound for the strength of the magnetic field in voids. Their claim, if true, 
is significant as it would represent the first evidence for magnetic fields in 
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such low-density regions and on such large scales. The basic idea is that an 
IGMF will deflect charged particles that arise in an electromagnetic cascade 
whose source is the very high energy 7-rays produced in an AGN. A non- 
observation of 7-ray secondaries coincident with VHE 7-rays, is then taken 
as evid ence of magnetic de flection along the path from the source to the ob- 
server. iNeronov and Vovkl dioio) quote a lower bound of i? > 3 x 10 G 
based on observations by the Fermi and the High Energy Ste reoscopic Sys- 
tem ( HESS) telescopes. A similar analysis was performed by lAleksic et al.l 
( 20101 ) using data from the MAGIC telescope though they are more cau- 
tious in presenting their conclusions and model dependencies. Note that 
iNeronov and Vovla ( 2010t ) assume continuous emission of gamma-rays for 10^ 
years, or longer which may not be realistic. If t his assumption is rel axed, then 
one obtains a more conservative lower bound ( Dermer et al.l[20TTh . 

The detection of the IGMF in filaments, if it exists, might be made 
with the next generation radio facilities. These facilities include the Square 
Kilometer Array (SKA), and upcoming SKA pathfinders, the Australian 
SKA Pathfinder (ASKAP) and the South African Karoo Array Telescope 
(MeerKAT), as well the Lo w Frequency Array (LOFAR) (see, e.g., papers in 
ICariUi and Rawlingsll2003 ). 

For discussions on ma gneti c fields in the ICM a nd cluster outskirts, see 
ISchekochihin et al.l (|2010[ ) and lBriiggen et al.l (|2010t) of this volume. 



3 Plasma Physics Mechanisms for Seed Fields 

3.1 Biermann battery 

The Biermann battery is a promising mechanism for the creation of astro- 
physical magnetic fields. The mechanism arises in an ionized plasma when- 
ever baroclinity exists, that is, when isodensity surfaces do not coincide with 
isobaric surfaces. This situation leads to an extra pressure gradient term in 
Ohm's law which drives currents. These currents, in turn, generate magnetic 
fields at a rate given by 

dt^e pI ' ^> 

wher e pe and Pe are the electron density and pressure, respectively (jBiermannl 
Il950|). Though t he mechanism was originally studied in the context of stars 



( BiermaimlllQSOf ) , it may arise during structure formation in cosmology when- 



ever the electr on pressure and density gradients are not aligned, as often oc 



curs i n shocks 

l2000LIXu et al.ll2008D . 



udritz and Sili3ll989llKulsrud et al ]ll997llDavies" and Widrowl 



We note that vorticity, w, is generated when the total pressure and density 
gradients are not aligned (baroclinity of flows - see Equation |23| below) . In 
an ionized plasma, a simple order of magnitude estimate yields 

B ^ "^u: ~ 3 X 10-21 ( . ^^——^ Gauss . (2) 

e \ km s ^ kpc ^ / 
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Since the present-day vorticity in the IGM is of order ^ a few km s^^kpc^'^ 
(see Section 15.11 below) , the seed fields from the Biermann battery will be 
rather small (~ 10~^° G). However, vorticity, and hence seed fields, were al- 
most certainly larger at early times. The following argument, based on dimen- 
sional analysis, suggests that for the objects which dominate the structure 
hierarchy at a given epoch the ratio of the vorticity to the Hubble parameter, 
H, is roughly constant and of order a few hundred. Once an object collapses, 
M ~ v^R/G where M, v, and R are the characteristic mass, velocity and 
size of the object. The spherical collapse model, which serves as a useful toy- 
model for structure formation in an expanding universe, predicts that the 
mean density within an object is a factor fc ~ 200 times the critical density, 
Pc- Moreover, objects tend to follow the mass-size relation M cx i?^. If we 
put all this together and use the fact that the critical density scales as the 
square of the Hubble parameter, we find that w ~ 300H (x {1 + zf^. Thus, 
we expect larger seed fields at early times. Similarly, the dynamical time for 
objects formed in the early universe is shorter than for objects today and 
hence the amplification of the magnetic fields by dynamo action will be more 
rapid. In short, the Biermann battery-dynamo amplification process may well 
operate at all stages of the structure-formation hierarchy. 



Various alternatives have also been considered. For example. [ Lazarian 



1992 ) considered a battery driven by electron diffusion. ISubramanian et al 



19941 ) considered the Biermann effect during the epoch of reionization when 



ionization fronts sweep through the medium generating currents and mag- 
netic fields. The fields are again fairly modest, but could be large enough to 
seed a dynamo which then amplifies them to an astrophysically interesting 
strength. 



3.2 Thermal fluctuations 

At finite temperatures, a plasma exhibits a finite (though possibly low) level 
of thermal fluctuations at all scales L = 2-K/k and frequencies w. Differences 
in the thermal motions between the different charges generate micro-currents 
and hence electromagnetic flelds which, on average, provide a noisy electro- 
magnetic background. In the early universe, the gas is both hot and dilute. 
For example, the gas temperature is T ~ eV both at the recombination epoch 
and after reionization but before the formation of LSS. Therefore, the cosmic 
gas constitute s a classical plasma to which the classical fluctuation theory 
(ISitenkol[l967h can be applied. 

According to this theory, the power spectral density of magnetic field 
fiuctuations b at wavenumber k and frequency lj is obtained from the spatial 
correlation function of the fiuctuating magnetic fields and is given by 

Here, the fiuctuating field is designated by lower case letters, i,j — 1,2,3, 
V? = k'^c^ /uP' is the index of refraction and e±{uj, k) is the complex transverse 
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dielectric response function of the isotropic plasma. The latter, in (nonrela- 
tivistic) thermal equilibrium, is given by 

where the sum is over negative (e) and positive (i) charges, uje.i, the respec- 
tive plasma frequencies, and Fq = dF^/dv, the derivative of the equilibrium 
distribution Fq which is assumed to be Maxwellian. The response function 
becomes 



e± = 1 



(5) 



Here, = {mi/mf,)9, and <l>{w) « 207^ + . . . is the Gordeyev integral of 
VD = y/3/2{uj/kcs) with = STe/rrie, the square of the sound speed, and 
6 = Te/Ti. Only at very large ion temperature T, ^ T^. do the ions contribute 
to the imaginary part. 

Seed fields for a turbulence dynamo arise from the zero frequency contri- 
bution, ci7 — >■ (index in the following equation) , which yields the magnetic 
power spectral density 



(&.&j-)0fc MO^eC^ e'/^kX, 10-'" VV 



(6) 



where = Te/nrieC^ is the normalized temperature and TefoV] is tem- 
perature in eV, Ae — c/wg ~ is the electron inertial length, N is 
the number density, and A'^[cm-3] the density in cm~'^. To obtain the last 
equality we have ignored the term fc^Ag which is negligible on galactic and 
extragalactic scales. The magnetic spectral energy density thus grows linearly 
with wavenumber fc, i.e. it decreases with the scale L. 

For comparison, a i? = 1 /xG field has the spectral energy density (-Bi^Q)ofc 
« 4.3 X 10^^* V^s'^/m. Were we to have integrated over all frequencies w, 
we would have obtained the equipartition between magnetic and electron 
thermal energy densities (&2)k — S/ioTe which holds in thermal equilibrium. 
Evidently, the thermal fluctuation level in Equation ^ is quite low. In order 
to obtain a B = 1/iG field, the spectral energy density must be amplified by 
the large factor of ~ 5 x lO'^^Ly^pcN^cm-'-^/VTcV- 



3.3 Filamentation instability 



Dynamos on galactic or cluster scales need seed fields on comparably large 
scales to get them started. Moreover, large-scale dynamos require turbulence 
and/or helicity which may not exist at adequate levels. Recently , the so-called 
filamentation (or Weibel) instability (IWeibellll959l lFriedlll959() has been dis- 
cussed as a generator of the IGMF (iGruzinov and Waxman E999I [g ruzmovl 
I2OOL .Medvedev and Loeb .199& ) and a possible alternative to dynamos on 
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galactic and extragalactic scales. The Weibel instability does require a strong 
departure of the plasma state from thermal equili brium, which can be pro- 
vided by fast beams on the plasni a background ()Achterberg and Wiersmal 
l2007l lFriedlll959i ISakai et aLll2004[ ) or temperature anisotropies. Th e latter 
can be caused by p ressure anisotropies fWeibel 'lOSQ', 'Treu mannll2010 h , shock 
waves (jjaroschek et al. 2005. Mshik awaLet^l , 2009), or thin current sheets. 
Turbulence, if it exists in the early universe, leads to thin curr ent sheets which 
can accelerate par ticles by second- order Fermi acceleration ( Jaroschek et al.l 
[2OO8, Jaroschek a nd Hoshindl2009l ). the local electric field which is generated 
in the current sheets, or by reconnection. Here, we discuss the possibility that 
the instability is driven by kinetic pressure/temperature anisotropics in the 
absence of magnetic field, as a mechanism to create seed fields. 

Pressure and temperature anisotropies are probably superior to the beam 
for driving the instability for the simple reason that beams are highly unsta- 
ble with respect to high-frequency plasma instabilities (Langmuir and Bune- 
man modes). If the beams have velocity Vf, and are current-compensated 
(i.e., the velocities of particles with different charges are the same) then the 
two-stream instability is excited with growth rate ^ts — V5{Nb/2Ny/'^u}e- 
Current-compensated beams arise in the symmetric counter-streaming con- 
figurations found in shocks. For the reasonable beam density Ni, ^ O.IN 
we find jts ~ 0.6 We ~ few kHz with a wavenumber is kts ~ <^e/Vf,. If the 
beams are not current compensated, which is the case in turbulence where 
a multitude of thin current sheets is generated, then the Buneman insta- 
bility is excited with growth rate 73 ~ 0.03 ~ 10^ Hz, if the current 
velocity Vc = Ve — Vi > Ve exceeds the electron thermal speed Ve- Other- 
wise, for Vc < Ve, the ion sound will be excited with growth rate roughly 
7ia ujpi ~ few 10^ Hz. All these instabilities grow very fast and readily de- 
plete the beams/streams, heating the electron plasma in the direction of flow 
and thus causing the pressure or temperature anisotropy of electrons with 
higher temperature along the flow direction, i.e. || to Vc. If this is assumed 
to be the (readily achieved) flnal microscopic state, then the temperature 
anisotropy A = Tgy/Te — 1 > becomes unstable with respect to the Weibel 
thermally anisotropic mode. The temperature anisotropy corresponds to a 
pressure anisotropy Pg — N[Te\ + {Te\\ — Te)VcYc/Vc] - One should note that 
only the electrons, because of their large mobility, are important; the ions are 
much less active on these microscopic time scales and thus do not contribute. 

In this thermally anisotropic case, the electrons obey a bi-Maxwellian 
equilibrium distribution function 



fe{vA 



(me/2^)3/2 



■ exp 



2T, 



2T, 



(7) 



and the Weibel instability takes over with the linear electromagnetic disper- 
sion relation = e± and the transverse dielectric response function 



1 



4o 



iA+i)[i + czm-^ 



(8) 



Here, Z{C,) is the plasma dispersion function, C, = uj/k±Ve, and Ve = \/2T'e/me 
is the electron thermal speed perpendicular to current. The Weibel instability 



8 




Fig. 1 The growth rate 7 — of the anisotropic-thermal Weibel instabihty nor- 
mahzed to the maximum growth, 7/7m, as a function of the normahzed wavenum- 
ber k/ko (see the text for the definition of ko). The vertical line indicates the 
position of the maximum growing wavenumber km/ko — l/v^- 

grows in the plane perpendicular to the direction of higher thermal velocity, 
which in our case has been assumed as the parallel direction. Hence, k = k^. 
The contribution of the resting ions has been retained for completeness; be- 
cause of the smallness of the ion plasma frequency LUi ^ We, being much less 
than the electron plasma frequency it plays no role in the instability but 
is important in the discussion of the thermal level. 

At zero real frequency, uj = 17, and the instability "f{kj_) > sets on for 
phase velocities uj/k± <^ Ve and wavenumbers k± < fco with growth rate 

vanishing at infinite wavelength k± ~ 0, where fco-^e = \/A. The growth rate 
maximizes at wavenumber k^Xe = fco-^e/'v/S = \J Aji (see Figure [TJ, where 
its value is 

7„ « 34y^iV[,„,-3]re[cV]^'/'(^ + 1) Hz. (10) 

A substantial growth can be achieved and stationary (purely growing) non- 
fluctuating, short-wavelength magnetic fields are produced. The maximally 
growing wavelength is indeed very short. For a plasma dens ity o f N « 0.01 
cm~'^, the wavelength is only of the order of A = (27rc/we)-y/3/A ^ 50\/3/A 
km. Even very weak anisotropics of A ~ 10"** give Am ^ 10* km. Galactic 
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or intergalactic scales L are far longer. For a realistic anisotropy of A <^ 1, 
at these large scales the growth rate scales with L as 

7iv ~ AaV^{c/L) ~ 10-"^v/7^/ikpc Hz (11) 

The above tells that mostly magnetic fields of small scales are excited 
and saturated by the instability, and those of galactic or intergalactic scales 
would not be substantial. Yet for ^4 ~ 0.1 and T ~ eV, the growth rate is 
7m/ ^ 10^^"* Hz on kpc-scales which is quite large. Magnetic fields of ^ 10^^^ 
G on kpc-scales would arise in ^ 10^ years or so provided the fields on smaller 
scales do not saturate the instability. 



4 Magnetic Fields after Recombination 

If magnetic fields were created before or during recombination they could 
have had a significant impact on the thermal and chemical evolution of gas 
during the dark ages, on t he formation of first stars and galaxies, and on the 
epoch of reionizat ion (e . g., Sethi and Subramanianl2005llTashiro and Sugivamal 
2006L ISethi et al.|[2008l ISchleicher et al.ll2008l |2009[ ) (see also lWidrow et al.l 



201d . of this volume). On the other hand, the first stars and galaxies may 
have been a source of magnetic fields. Our discussion should be viewed with 
a critical eye since the exis tence, strength and scale of e arly universe fields 
is highly uncertain (see e.g. iGrasso and Rubinstein|[200ll) as is the efficiency 
of dynamo action at early times. 



4.1 Magnetic fields during the dark ages 

The dark ages designate the period between recombination and reioniza- 
tion, during which no astronomical objects exist. The only radiation pro- 
duced during this epoch is the 21 cm line of neutral hydrogen. The pri- 
mordial gas consists of about 75% hydrogen and 25% helium, as well as 
trac eable amounts of lithiu m and beryllium (both of the order 10~^) (see, 
g., Kolb and Turiieil 19901). The ionization degree free zes out at about 10" 



(|Peebleslll968llZeldovich et al.lll969llSeager et al.lll999D . Though this ioniza- 



tion degree may seem rather low, it is significantly higher than the typical 
ioniz ation degree in present-day molecular clouds (jMvers and Khersonskvl 
Il995[ ). More to the point, it is high enough to prevent non-ideal magneto- 
hydrodynamic (MHD) processes such as ambipolar or Ohmic diffusion from 
dissipating the magnetic fields. 

The more relevant question is whether magnetic fields can be sustained 
during gravitational collapse when both the density and recombination rate 
increase and hence the ionization degree decreases. Gravitational collapse 
occurs on length scales where the thermal pressure can no longer balance the 
gravitational force. The critical length scale is given by the Jeans length 

2\ 1/2 
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wh ere c., is the thermal so und speed and p the mass den sity. Previous studies 
by iMaki and Susal (|2004l) and iGlover and Sayinl (|2009[ ) followed the chemi- 



cal evolution during gravitational collapse. Their one-zone model takes into 
account the density evolution in the central core, where the collapse takes 
place on the free-fall timescale i// ~ Xj \JGp. As is well-known from analytic 
solutions and numerical simulations, the density i n the central core is roughh 



homogeneo us on the scale of the Jeans len gth Aj 
lAbel et allbool iBromm and Loeb"2003', 'Yoshk 



(iLarsoii 



Yoshida et al. 



atrai core is rqugniy 
1969l . lPensto3l 19691 



2008). 



The calculations of ,Maki and Susa ( 2004) and iGlover and SavinI (|2009f l 
showed that the ionization degree evolves roughly as p~^l'^ until a density 
of ^ 10^ cm~^ is reached. At this stage, the recombination rate becomes so 
high that the proton abundance becomes negligible. However, the presence 
of lithium and the low recombination rate of Li"*" stabilizes the ionization 
degree at a level of ~ 10~^°, and it stays roughly constant even at higher 
densities. 

In order to assess the implicat ions of this result. ISchleicher et al.l (|2009f ) 
combine the chemical network of [gi over and SavinI (|2009f ) with a state-of- 
the-art model for am bipolar and Ohmic dif fusions, based on the work of 
iPinto et all (|2008t) and lPinto and Gallil (|2008[) . In their multi- fluid approach, 
the ambipolar diffusion rate is given as 

iAD = ^|(VxB)xB/S^ (13) 
47r 

where the B denotes the magnetic field. The ambipolar diffusivity t^ad is 
given as 

'^Ai -II'/Ai,z- (14) 

i 

with 77AD.i denoting the ambipolar diffusivity due to an ionized species i. 
Similarly the Ohmic diffusion rate is given as 

iohm = 7^ |V X B|^ (15) 
47r 

with the Ohmic resistivity t]o given as 

^o^Y.<^- (16) 

i 

The ambipolar resistivities ?7AD,i and the Ohmic resistivities r/o,i due to ion- 
ized species i are calculate d from the m omentum-transfer coefficients as de- 
scribed in the appendix of iPinto et al ] (12008) . 

To guide the physical intuition, we briefly summarize their results for the 
case of two charged fluids, e.g. electrons and positive ions, but note that the 
more detailed multi-fluid approach has been adopted in the calculations. The 
key quantities that regulate non-ideal MHD effects are the Hall parameters 
/3sm, dcflned as 

/3sn=f^)^^^. (17) 

Here, Qs denotes the charge of the species s, its mass, m„ the mass of 
the neutrals, p„ the neutral mass density, c the speed of light and {av)sn the 
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momentum-tra nsfer coefficients betwe en the ions and the neutrals (for the 
ex pressions, see 
bv lPinto et all I 
as 



Pinto and Gallil ( 20081 ) and the discussion below). As shown 



2008f ). the ambipolar and Ohmic resistivities are then given 



VAD = 



Vo = 



\ A-Ka J (3. 



1/3. 



(18) 
(19) 



where f3~n denote the Hall parameters for the positively and negatively 
charged species, and a is given as 



and 



(20) 



the mass 



In the latter expression, denotes the positive charge 
density of the positively charged species. 

For colli sions in volving protons and electrons, the momentum-transfer co- 
efficients of |Pinto"an d Galli (200^ are adopted. For collisions involving Li+ , 
the momentum-transfer coe fficients using the polar ization approximation are 
calculated, as described by 'Schleicher et al." (2009'). This is indeed justified 
for collisions with helium (jCas sidv and Elford 1985) and H2 (Dickinson et al] 
Il982l iRoeggen et al.ll2002[ ). For collisions with atomic hydrogen, no detailed 
theoretical or experimental measurements are currently available. This pro- 
cess is of minor importance, however, as Li"*" becomes the dominant ionized 
species at densities when molecular hydrogen is the dominant neutral species. 

Equations (|13p and (|15p show that the dissipation rates depend on the 
magnetic field strength itself, and we note that also the ambipo lar resis- 
tivities increase with increasing field strength (jPinto et al ] |2008D . During 
gravitational c ollapse, an initi ally weak field, for instance due to the Bier- 
mann battery (|Xu et al.ll2008l) . may be ampl ified both due to grav itational 
compression as well as turbulence dynamo ( Schleicher et al.l l2010[ ) . If the 
amplification due to turbulence dynamo is efficient, s aturation may occur 
at equipartition or a somewhat lower level 
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.Schekochihin et al.l 12004 . 
I2009D . As ambipolar and 
Ohmic diffusions are stronger for stronger fields, we will initially focus on 
equipartition magnetic fields and show that magnetic fields can be sustained. 
The same is then true also for weaker magnetic fields. 

To have a definite model, we assume Kolmogorov turbulence, and thus 
that the equipartition field strength scales as l-^^^. This is a good approx- 
imation for subsonic turbulence where compressional effects are of minor 
importance, and it is this type of tur bulence that is expected in the first 
star- forming systems ( Abel et"al]l2002[ ). The Ohmic and ambipolar diffusion 
scales are calculated as a function of the central core density and compared 
them to the Jeans length (see Figure [5]). For clarity, the Ohmic / ambipo- 
lar diffusion scales are defined as the length scales on which the Ohmic or 
ambipolar diffusion time equals the eddy-turnover time. The calculation as- 
sumes that the turbulent velocity and the magnetic field follow a typical 
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Fig. 2 Diffusion length scales during gravitational collapse. The black solid line 
denotes the Jeans length as a function of central core density. The ambipolar diffu- 
sion scale for the equipartition magnetic field is given by the red dashed line, while 
the red dashed-dotted line corresponds to the scale for a weaker field (the magnetic 
energy equals 1/60 of the kinetic energy). The blue dotted line gives the Ohmic 
diffusion scale, which is independent of the magnetic field strength. In all cases, the 
dissipation scales are smaller than the Jeans length. 



scaling law as expected for Kolmogorov turbulence. During collapse, the in- 
creased densities lead to a more efficient coupling between the ions and the 
neutrals, which largely compensates for the decreasing ionization degree. At 
a density of 10^ cm~^, the ambipolar diffusion scale increases due to a drop 
of the ionization degree at these high densities, but still stays an order of 
magnitude below the Jeans length. The Ohmic diffusion scale always stays 
orders of magnitude below the Jeans length. The figure also shows that when 
the magnetic field strength is below the equipartition value, the ambipolar 
diffusion scale is even small. The Ohmic dissipation scale, however, stays on 
the same low level, as the Ohmic dissipation rate scales roughly with the 
magnetic energy. 

Numerical calculations suggest that the ionization degree may even be 
higher than expected from one-z one models, in particular due to the p resence 
of shocks of Mach number ~ 1 (jClark et al.ll201ll iTurk et ahlboogf ). So we 



safely conclude that during the dark ages, the magnetic fields up to the 
equipartition strength can be sustained even at high densities in collapsed 
regions. 
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Fig. 3 The magnetic field strength after a free-fall time. The horizontal axis denotes 
the length scale of the turbulent eddies. We assume a power-law scaling relation for 
the rms turbulence on a given length scale, and explore the cases of Kolmogorov 
and Burgers-type turbulence (Schleicher et al. 2010). In this model, the turbulent 
injection scale is 100 pc. For every length scale, we show both the maximum field 
strength that could be obtained from the small-scale dynamo, as well as our esti- 
mate for the actual field strength that might be produced within a free-fall time. 
The model suggests the production of strong magnetic fields on small scales. 



4.2 Magnetic field evolution in first star-forming halos 

We now turn to the amplification of magnetic fields during the formation of 
the first stars. The halos that harbor the first stars at redshifts of z ~ 20 — 30 
typically have a total mass of 10^ Mq, consist of primordial gas with a 
temperature of ^ 3000 K, and have a spatial extent of ^ 100 pc. Turbu- 
lence is expected to be present, but on a subsonic level ()Abel et al.ll2002l 
iBromm and Loebll2003l . f^shida et al.. .20081 Such turbulence may create 
strong magnetic fields in the first galaxies and around them (iRvu et al.ll2008L 



lArshakian et al.ll2009l Tde Souza and Opheij[20Toi ISchleicher et al.ll201Clir 



Consider a simple model in which magnetic fields are amplified and satu- 
rate at a level corresponding to an energy 1/60 of the equipartition energy on 
an eddy turnover time ioddy We assume that turbulence is injected during 
the virialization phase on the length scale of the system (i.e., ^ 100 pc) with 
the velocity of the order of the sound speed. Both Kolmogorov and Burgers 
turbulences are considered; for Kolmogorov turbulence the velocity depends 
on the length scale as l^/^ , while for Burgers turbulence the scaling accords 

to ;i/2. 

For magnetic fields to be relevant during the formation of the first stars, 
they must be amplified within the free-fall time scale i// ^ Xj \JGp. The 
number of eddy turnovers is thus given by the ratio tff /teddy Here, the 
eddy turnover time teddy = l/v has a characteristic dependence on length 
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scale for Kolmogorov and Burgers turbulences. Then, the expected magnetic 
field strength is given in Figure [3] as a function of the length scale of tur- 
bulence. The preliminary calculation shows that magnetic fields of the order 
10^^ G can be reached in the first star-forming halos. For a comprehensive 
understanding of the strength and structure of the magnetic fields, the de- 
tailed implementation of turbulence dynamo will be required. Recent numer- 
ical MHD simulations suggest that the small-scale dynamo indeed operates 
during gravitational co llapse, pro ducing mag netic fields durin g the formation 
of the first structures (|Sur et ah 2010, Fcdcrr ath eraDl201ll) . 

Additional amplification may occur after the formation of protostellar 
disks. In particular, a large-scale dynamo and the magneto-rotational insta- 
bility may further enhance the magnetic fields in the first star-forming halos 
(jPudritz and^iigil989llTan and Blac kman"2004'.'Silk and Langeil l2n06V For 
simplified field geometries. iMachida et ah, (20Q(i, iOOSi) found that magnetic 
fields lead to the formation of jets and help to suppress fragmentation in 
protostellar disks. IPVomang et al.l ( 2004 ) studied self-gravitating magnetized 
disks and found that the interaction of turbulence created by the magneto- 
rotational instability may excite additional modes for the gravitational insta- 
bility, and that the interaction of these modes reduces the accretion rate. This 
work suggests that magnetic fields may be amplified rapidly in the first star- 
forming halos and may become dynamically relevant though more realistic 
studies are necessary. 



4.3 Seed fields from astrophysical processes 

The first stars probably possessed strong magnetic fields and therefore may 
have provided seed fields for dynamos in galaxies and in the IGM. If the 
stars subsequently explode as supernovae or lose a significant amount of 
mass through stellar winds, the fields ejected along with mass will find their 
way into the interstellar medium (ISM) and spr ead beyond galaxie s into the 
IGM through galactic winds. Simple estimates bv lSvrovatskii ( 1970t ) illustrate 



the viability of the process. If there have been some 10 supernovae over 
the lifetime of galaxies, each of which spreads material through a (10 pc)^ 
volume. Using values for the field strength typical of the Crab nebula, one 
therefore expects galaxies to be filled by 10 pc regions with fields of strength 
^ 3 /iG. Assuming the same L~^l'^ scaling, one finds fields of strength of 
~ 10^^^ G on 10 kpc scales. This value is si gnificantly larg e r tha n those 
from the processes described in Section |3] and in lWidrow et al. (l2010l) of this 



volum e, although the filled volume is rather small. Recently, Donnert et al.l 



( 2009f ) suggested that galactic outflows during the starburst phase of galactic 



evolution can deposit a substantial amount of magnetic fields in the IGM. 

Seed fields can be produced at cosmologic al shocks which were induced 
durin g the format ion of the LSS of the univers e (iRvu et al. 12003 . IPfrommer et al.l 
l20nfiL iKang et al.ll2007l ISkiUman et al]l2008l IVazza et al.l l2009l (see Section 



15.11 for further description of cosmological shocks) . Cosmological shocks are 
collisionless like shocks in other astrophysical environments, where CRs are 
accelerated at the same time as the gas is heated. During the process of ac- 
celeration, it was shown that the upstream magnetic field can be amplified 
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nonlinearly by non-resonant growing modes ( BeUll200"3 ). Then, the magnetic 
field can have the energy up to 

£b ^ -z £CRs- (21) 

2 c 



With Uni/c ~ 10 and Mach number a few in cosmological shocks ( Rvu et all 



l2003llKang et al.ll2007i) . we get es ~ ecRs ~ £therm- In addition, the Weibel 
instability described in Section 13.31 can operate and produce magnetic fields 



up to the level of e_B ~ 10 Csh (jSchhckeiser and Shuklall2nn3llMedvedev et alJ 



|2006[ ). Here, e^h is the energy density of upstream flow, and Sb, sgrs, and 
Ethcrm a^c the energy densities of downstream magnetic fields, CRs, and gas 
random motion, respectively. These processes can potentially produce strong 
fields around cosmological shocks, although the volume filling is small and 

the coherence length of the resulting fields is expected to be mi croscopic. 

Ot her mechanisms have been proposed. For instance, recentlv. lMiniati and Belli 
()201l[ ) suggested that the return current which is induced by cosmic-rays pro- 
duced by early supernovae can deposit seed fields into the IGM. 



5 Magnetic Fields in the IGM 

With the processes described in Sections [3] and |4] and also in IWidrow et al.l 
(l20inl) of this volume, there is no shortage of mechanisms to generate seed 
fields for the IGMF. Those fields are expected to be amplified by the turbulent 
fiow motions which were induced during the formation of the LSS of the 
universe. The turbulence dynamo not only increases the strength of magnetic 
fields, but also produces the magnetic fields of large scales, up to the energy 
injection scale, through the inverse cascade. The beauty of the turbulence 
dynamo in LSS is that it erases the memory of weak seed fields and produces 
the IGMF, independent of the origin of seed fields. 



5.1 Turbulence in the LSS of the universe 



Signat ures of turbulence have been observed in the ICM. For example. lSchuecker et al.l 
(|2004D analyzed gas pressure maps that were constructed from XMM-Newton 
X-ray data. They claimed that in the Coma cluster, which appears to be in 
a post-merger state, pressure fiuctuations are consistent with Kolmogorov 
turbulence. The turbulence is likely subsonic but with an energy that is 
at least 10 % of the thermal energy, i.e., eturb > O.lethorm- The results 
agree with predictions from numerical simulations, namely that the flows 
i n cluster scales have a power spectrum expected for Kolmogorov turbulence 
( Kulsrud et al 1 119971 iKulsrud and Zweibeil2008l^. and even in relaxed clus- 
ters the flow motions have etin ~ O.lethor,^ (^Na gai et al.l [2007). Tu rbulence 
in the ICM was also studied in RM maps of a fe w clusters ( Vogt and EnfilinI 
l2005llGuidetti et al.ll2008L iBonafede et al.ll2010t ). 

Recently, Rvu et al. (|2008[ ) proposed a scenario in which vorticity is gener- 
ated directly or indirectly at cosmological shocks and turbulence in the IGM 
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Fig. 4 Two-dimensional slice of (85 h ^Mpc)^ showing shock locations with color- 
coded shock speed (left panel) and the magnitude of vorticity (right panel) at 2: = 0. 
The vorticity is given in units of 10"'' tagl, where tage is the age of the universe. 



is induced via the cascade of the vorticity. Here, we provide an estimate esti- 
mate of the turbulence seen in simulations fo r the formation of t he LS S. The 
results are based on a simulation described in lCen and Ostrikeil (120061^ which 



includes the radiative proces ses of heating/coo ling and feedback from galac- 
tic superwind. The work of iRvu et al.l (200l) utilized a simulation, where 



only the gravitational and gas dynamical processes are included. 

In the simulation of ICen and Ostrikeil (j2006t) . the WMAPl-normalized 
ylCDM cosmology was employed with the following parameters: fib — 0.048, 
n„i = 0.31, Ha = 0.69, h = i?o/(100 km/s/Mpc) = 0.69, erg = 0.89, and 
n = 0.97. A cubic box of comoving size 85 Mpc was simulated us- 
ing 1024^ grid zones for gas and gravity and 512'^ particles for dark mat- 
ter. It allows a uniform spatial resolution of Al = 83 h~-^ kpc. For de- 
tailed descriptions of input physics ingredients such as non-equilibrium ion- 
ization/cool ing, photoi o nizati on/h eating, star format i on, a nd feedback pro- 
cesses, refer ICen et al ] (|2003D and [C en and Ostrikeil (12006). Feedback from 



star formation was treated in three forms: ionizing UV photons, galactic 
superwinds, and metal enrichment. Galactic superwinds were meant to rep- 
resent cumulative supernova explosions, and modeled as outflows of several 
hundred km s"^. The input of galactic superwind energy for a given amount 
of star formation was determined by matching the outflow velocities com- 
puted for star-burst galaxies in the simulation with those observed in the 
real world. The simulat ions were perfor med using a PM/Eulerian hydrody- 
namic cosmology code ( Rvu et al.lll993f ). 



In the IGM, vorticity, a; = V x v, can be generated directly at curved 
shocks and also by the baroclinity of flows. For uniform upstream flow, the 
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vorticity produced behind curved shock surface is 



P1P2 



KVi X h, (22) 



where pi and p2 are the upstream and downstream gas densities, respectively, 
Ui is the upstream flow velocity in the shock rest frame, K is the curvature 
tensor of the shock surface, and n is the unit vector normal to the surface. If 
isopycnic surfaces do not coincide with isobaric surfaces, vorticity is generated 
with the rate given by 

t^bc = ^ Vp X Vp. (23) 

P 

Shock waves are ubiquitous in the IGM, as in other astrophysical envi- 
ronments. The spatial distribution and properties of cosmological shocks in 
the LSS of the universe have been studied quantitatively using simula tions 
for the formation of LSS (iRvu et al.ll2nn3llPfrommer et all 2006. Ka ng eTall 



l20n7L ISkillman et al. I I2OO8I IVazza et al.ll2009l^ ." In the cold dark matter uni- 
verse with cosmological constant (vlCDM), shocks with Mach number up to 
a few hundreds and speed up to a couple of thousand km s^^ are present 
at the present universe (z = 0). In the left panel of Figure HI the spatial 
distribution of cosmological shocks is shown. Numerous shocks are found. 
External shocks exist around sheets, filaments, and knots of mass distribu- 
tion, which form when the gas in void regions accretes onto them. Within 
those nonlinear structures, internal shocks exist, which form by infall of pre- 
viously shocked gas to filaments and knots, and during subclump mergers, as 
well as by chaotic fiow motions. Due to the low temperature of the accreting 
gas, the Mach number of external shocks is high, extending up to Af ~ a few 
xlOO, while internal shocks have mostly low Mach number of M ~ a few. 
The mean distance between shock surfaces is ~ 3 h~^Mpc when averaged 
over all the universe, or ^ 1 h~^M.pc inside nonlinear structures. Internal 
shocks of M ~ 2 — 4 formed with hot and high-density gas are responsible 
for most of shock dissipation into heat and CRs. It was shown th at the shock 
dissipation can count most of the gas thermal energy in the IGM ( Kang et al.l 
[2005h . 

In the right panel of Figure |H the distribution of vorticity is shown. 
It closely matches that of shocks, suggesting that a substantial portion of 
the vorticity, if not al l, has been generated at the shocks. As a matter of 
fact, as was noted in iRvu et aD ( 2008( ). the vorticity in the IGM can be 



accounted with that generated either directly at curved cosmological shocks 
or by the baroclinity of fiows. The contributions from the two processes are 
comparable. The baroclinity resulted from the entropy variation induced at 
shocks. So all the vorticity generation also can be attributed to the presence 
of cosmological shocks. 

For quantification of vorticity, the flow velocity is decomposed into 

V = Vdiv+Vcurl+Vunif, (24) 

where the divergence and curl components are deflned as V • vjiv = V • v 
and V X Vcuri = V x v, respectively. That is, Vdiv is associated to compres- 
sional motions, while Vcuri to incompressible shear motions. Here Vunif is the 
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log(ic/27T) (l/85h-iMpc) 

Fig. 5 Power spectra, J Pydk = {^{l/2)v'^^, for the gas velocity and its curl and 
divergence components at 2 = 0. Two straight lines of slopes —5/3 and —1 are also 
drawn for comparison. 



component uniform across the computational box, whose magnitude is much 
smaller than the other two components. The decomposition is calculated 
exactly in Fourier space. We note with the above decomposition, locally 

Vdiv • Vcurl 7^0 so \ (^^div + ^'curl) ■ (25) 

However, globally 

/ Vdiv • Vcurlrf^X = so / ^W^d^X = \ \ (Wdiv + ^'curl) '^^^ (26) 
^box v/box ^ v/box ^ 

The power spectra for the gas velocity and its curl and divergence com- 
ponents at the present universe are shown in Figure [S] At long wavelengths, 
the amplitude of perturbations are small, so that linear theory applies. That 
is, Pcn^xik) — >■ as fc — )• 0, while Pdiv(fc) follows the analytic theory expecta- 
tion, Pdiv(fc) ^ . For wavelengths smaller than a few Mpc, nonlinearities 
dominate, and we see Pcuri(fc) ^ -Pdiv(fc)- ^curi(fc) peaks at ~ 5 Mpc, 
and for fc somewhat larger than the peak wavenumber, the spectrum follows 
a power law of fc~^/'^, the Kolmogorov spectrum. PcurK^) has most power at 
~ 2 — Zh^^ Mpc, that indicates the typical scale of nonlinear structures in 
the simulation. 

The vorticity in the IGM as a function of gas density and temperature 
is shown in the left panel of Figure [S] The figure exhibits a clear trend 
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Fig. 6 Distribution of the magnitude of the vorticity in the IGM (left panel) and 
the strength of the IGMF (right panel) in the density-temperature plane at z = 0. 
The vorticity is given in units of f^ge- 

that the vorticity is larger in hotter regions. It is because hotter regions 
are occupied with the gas that has gone through shocks of larger speed. The 
vorticity generated at shocks of larger speed should be on average larger. Also, 
there are regions of high density (pgas/ (Pgas) ^ 100) and warm temperature 
(T ^ 10'* K), where the vorticity is large. These regions contain the gas that 
was heated to high temperature and subsequently cooled down. 

In the left panel of Figure [71 the root-meas-square (rms) of the vorticity 
at a few redshifts for the gas in the four phases of the IGM is shown. Note 
that the vorticity shown was computed over the same comoving scales, and 
normalized with the age of the universe at given z. The vorticity increases 
as the universe evolves and the LSS of the universe develops. But over the 
period of the time presented in the figure, the vorticity, especially in the hot 
IGM and WHIM, has increased just by a factor of a few. 

Figures E] and [7] indicate that at the present epoch, Wrmsiagc ~ a few xlO 
inside and around clusters/groups with T ^ 10^ K and - 10 in filaments 
which is filled mostly with the WHIM. On the other hand, uj^mstage is on 
the order of unity in sheetlike structures and even smaller in voids. Here, 
tagc is the present age of the universe, so Wi-ms^agc represents the number of 
eddy turnovers in the age of the universe. It takes a few turnover times for 
vorticity to decay and develop into turbulence. So it is likely that the flows in 
clusters/groups and filaments is in a turbulent state, while turbulence have 
not significantly developed in sheetlike structures and voids. 
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Fig. 7 The time evolution of the rms of the vorticity for four temperature phases 
of the IGM (left panel) and the time evolution of the averaged strengths of the 
IGMF for the hot IGM and WHIM (right panel) as a function of redshift z. The 
vorticity is given in units of tags (a) where tage(2) is the age of the universe at z. 
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Fig. 8 Turbulence to thermal energy ratio as a function of temperature at « = 0. 
The values shown are volume-averaged and mass-averaged over temperature bins. 
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To estimate the energy associated with the turbulence induced in the 
IGM via the cascade of the vorticity, we assume that the energy of vortical 
motions, (l/2)/5gas'^curi' is transferred to that of turbulent motions. Then, 
we can regard it as the turbulence energy, eturb- Figure [5] shows the ratio 
of the turbulence to thermal energies in clusters/groups and filaments as 
a function of temperature at the present universe. In clusters/groups with 
T <; 10^ K, Eturb < £thorm- Particularly the mass averaged value shows that 
^turb/^thorm would be ~ 0.1 — 0.3 in cluste r cores. The predicted £turb/£thcrm 
is in a good agreement with observation (jSchuecker et al ][200l. Note that 
Mturb = Wcuri/cs = v^TS (eturb/ethcrm) , whcrc Cs IS the sound speed. 
Therefore, overall turbulence is subsonic in clusters/groups, whereas it is 
transonic or mildly supersonic in filaments. 



5.2 Amplification of the IGMF by the turbulence in the IGM 

In principle, if simulations for the formation of LSS in the universe includes 
magnetic fields, that is, if they are MHD, the amplification of the seed mag- 
netic fields by the turbulent motions in the IGM should be able to be fol- 
lowed. But, in reality, simulations with the current capacity of computing 
power have too low a resolution to reproduce the full development of turbu- 
lence inside nonlinear structures. Also the numerical resistivity is larger than 
the physical resistivity by great many order of magnitude. As a result, the 
growth of magnetic fields is saturated before dynamo action become fully op- 
erative, and the amplificati on of magnetic fields can not be followed correctly. 



as was already pointed in (jKulsrud et al.i 19971 



In order to reproduce the growth of magnetic fields by dynamo action, a 
separate three-dimensional simulation of MHD turbulence in a controlled 
box was performed; incompressible, driven turbulence with initially very 
weak or zero magnetic fields was simulated using a pseudospectral code 
( Cho and Vishniacl2000t) . Hyperviscosity and hyperresistivity with the Prandtl 



number of unity were used. The advantage of performing incompressible sim- 
ulations using a pseudospectral code is that the intrinsic numerical viscosity 
and resistivity are virtually zero. And by using hyperviscosity and hyperre- 
sistivity, the inertial range can be maximized. We point that it would take 
much higher resolution to achieve the same growth rate in simulations of 
compressible MHD turbulence. 

The top-left panel of Figure IH] shows the time evolution of kinetic and mag- 
netic energies for four different simulations: 256H8-BolO^'^, 256H3-BolO~^, 
256H8-BoO, and 256H8-BolO-s. Simulations are denoted with 256Y-BoZ, 
where 256 refers to the number of grid points in each spatial direction, Y 
refers to hyperdissipation (H) and its order, and Z refers to the strength 
of the external magnetic fields. The turbulence was driven at the scale of 
La ~ (l/2)Lbox where Lbox is the computational box size, and the driving 
strength was set so that the total energy is i^tot = ^-kin + Eb ~ 1 at satura- 
tion. The time is given in units of the eddy turnover time that is defined as 
the inverse of vo rticity at driving scale, icddy = l/wdriving, at saturation. See 
[Cho et al.1 (|2009t ) for details of simulations. 
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Fig. 9 Top-left panel: The time evolution of V'^ and for four different runs. Hero, 
the kinetic energy and magnetic energy densities arc V'^ /2 and /2, respectively. 
Top-right panel: Time evolution of peak scales of magnetic field spectrum, L^^k), 
and LE(k)/k- Bottom-right panel: Time evolution of the energy equipartition 



Lk 



kE(k); 

scale, icq. Bottom-left panel: Time evolution of the integral scale, Lint, and the 
curvature scale, Lcurv ^'s are the scales normalized with the energy injection scale 
Lq. See the text for details. 



The amplification of magnetic fields by turbulence dynamo is shown in the 

figure. It does not sensitively depend on the initial magnetic field strength 
once it is sufficiently weak as well as details of simulations including the 
dissipation prescription. The evolution of magnetic fields goes through three 
stages: the initially exponential growth when the back reaction of magnetic 
fields is negligible, then the linear growth when the back reaction starts to 
operate, and the final sat\.iration. By fitting the evolution, we model the 
growth and saturation of magnetic field energy as 

r 0.04 X exp [(t/teddy - 4)/0.36] for t/teddy < 4 
Eb = < (0.36/41) X {t/Uddy - 4) + 0.04 for 4 < t/t^My < 45 (27) 

[ 0.4 for t/teddy > 45 

Along with the amplification of magnetic field strength, the magnetic 
fields become coherent through the inverse cascade. For the quantification, 
different characteristic lengths of magnetic fields can be defined in MHD 



23 



turbulence: the peak scale of the spectrum of magnetic fields, L_E(fc)) the scale 
containing the largest energy of magnetic fields, L]~E(k)i the peak scale of the 
spectrum of projected magnetic fields, L^f^f^yf^, the energy equipartition scale, 
Loq(= 27r/fcoq), defined as 



/ E^{k) dk= Eb{k) dk, (28) 

Jk^q Jo 

the integral scale, Lint, defined as 

jE,{k)/k dk 
= jE,ik)dk ' 

and the curvature scale, Lcurv, defined as a typical radius of curvature of field 
lines. 

The rest of Figure [S] show the time evolution of characteristic lengths of 
magnetic fields. At saturation, the peak of magnetic field spectrum, LE(k), 
occurs at ~ Lo/2, where Lq is the energy injection scale, while the most 
energy containing scale, L^^^ik) is ^ Lg/b. During the stage of magnetic field 
amplification, the energy equipartition scale, Leq, shows a power-law increase 
of ~ t^'^, while the integral scale, Lint, and the curvature scale, Lcurv, show 
a linear increase. The equipartition, integral, and curvature scales satu rate 
at ^ Lq, ^ O.SLq, and ^ O.ISLq, respectively. See ICho and Rvul ( 2009t) for 



further details of characteristic lengths in MHD turbulence with very weak 
or zero mean magnetic fields. 

The results of the incompressible MHD turbulence simulation were con- 
voluted to the data of the stimulation for the formation of the LSS of the 
universe to get the IGMF. For the estimation of the magnetic field strength, 
it was assumed that a fraction of the turbulence energy is converted into the 
magnetic energy. The fraction was expressed as a function of the number of 
local eddy turns over the age of the universe, so 

£B = Ecurl • X iago)- (30) 

For the fraction (/) (a; x tago), the fitting formula in Equation ^7} was used. 
Then, the magnetic field strength was calculated as (STres)^/^. 

The resulting magnetic field strength as a function of gas density and 
temperature is shown in the right panel of Figure IHl On average, the IGMF 
is predicted to be stronger in hotter and denser regions. It is because with the 
turbulence energy is larger in hotter and denser regions. Also the conversion 
factor, (f>, is larger in hotter and denser regions. 

In the right panel of Figure [71 the volume- averaged and rms values of the 
magnetic field strength at a few redshifts for the gas in the hot IGM and 
WHIM are shown. Our scenario predicts that the magnetic field strength 
would be (B) <; 1 jiG inside clusters/groups, ^0.1 /iG around clusters/groups, 
and ~ 10 nG in filaments at the present universe. The magnetic fields should 
be much weaker in sheetlike structures and voids. But as noted in Section 
15.11 turbulence is not fully developed there. So our estimation of the IGMF 
in those regions should not be applicable. 
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In each temperature range, the magnetic fields were stronger in the past at 
z ~ a few. Our calculation indicates that the vorticity has not changed much 
since z ~ a few ( see the left panel of F igure [71), nor the vortical component 
of flow velocity (iRvu and Kand [20081 ) . So the stronger magnetic fields in 
the past should be due to the higher density. We point that clusters were 
virialized around 2 ~ 1 or so, so the density within the virial radius of ~ 1 
Mpc has not changed much since then. But our estimation of the magnetic 
field strength for the hot IGM with T > 10^ K also includes the gas in cluster 
outskirt, which exte nds far beyond the virial radius, up to several Mpc or 
even larger (see, e.g. lRvu et al.ll20d3l) : the density of the gas there was higher 
in the past. We also note that the magnetic fields, when averaged all over 
the computational box, were weaker in the past, because the fraction of the 
strong field regions was smaller. 

Based on a kineti c theory, assuming th e K olmogorov spectrum for tur - 
bulence flow motions. iKulsrud et afl ()1997[ ) and lKulsrud and Zweibel (|2008[ ) 
estimated that the magnetic field strength in clusters would be a few fiG. 
Our result agrees with the previous work. It also matches well with the ob- 
served strength of magnetic fields in the ICM, which is discussed in Section 
[TJ On the other hand, our prediction of i? ~ 10 nG in filaments is within but 
lower than the upper Umit of ~ 0.1/xG whi ch is imposed from RMs outside 
of clusters (|Rvu et aLlllOOl IXu et al.ll2006[ ). 

The characteristic lengths of the IGMF in our scenario can be conjectured. 
In clusters, turbulence is near the saturation stage with t/ teddy ^ 30 (see 
Figures [7] and O . Then, for instance, we may estimate that the peak scale 
of magnetic field spectrum and the integral scale would L^^^-j/Lo ~ 0.4 
and Lint/Lo ^ 0.2. If we take the energy injection scale Lq ^ 100 kpc, 
which is approximately the scale height of cluster core, LE{k) ~ 40 kpc and 
Lint ~ 20 kpc, respectively. In filaments, on the other hand, with t/t^ddy ~ 10, 
turbulence is expected to be still in the linear growth stage, and LE(^f^j/Lo 
1/10 and Lint/Lo ~ 1/15. We may take the energy injection scale Lq ^ 5 
Mpc, which is the typical thickness of filaments. Then, LE(k) ~ 0.5 Mpc and 
Lint ~ 0.3 Mpc, respectively. 



5.3 Contribution from AGNs 

There is a possibility that the IGMF is further strengthened by the magnetic 
fields ejected through jets from back holes in AGNs (see, e.g.. lKronberg et al.l 
l200l[ ). Strong magnetic fields almost certainly arise in accretion disks sur- 
rounding black holes. These fields may find their way into the IGM via mag- 
netically dominated jets. The pote ntial fi eld strength due to this process can 
be estimated as follows (see, e.g.. lHovl3 il969). The rotational energy asso- 
ciated with the central compact objects of mass M which power AGNs can 
be parametrized as fMc^ where / < 1. If we assume equipartition between 
rotational and magnetic energies within a central volume Vc, we find 




(31) 
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If this field then expands adiabatically to fill a volume V in the IGM, one 
finds B - Bc{Vc/Vf^^. Considering the values of M = lO^M©, / = 0.1, 
and y ~ (1 Mpc)'^, we get B ^ 10 nG. Note this field strength is comparable 
to that estimated in filaments with turbulence dynamo. So this process may 
enhance the strength of the IGMF by a factor of two or so. 

Although there are some RM observations that indicate strong magnetic 
fields in some of jets, however, it is not yet clear whether all jets are mag- 
netically dominated. Also the details of the population of AGN jets and the 
volume filling fraction of magnetic fields in the universe by this process need 
to be further worked out. 



6 Astrophysical Implications of the IGMF 

The existence of the IGMF, especially in filaments, can have noticeable im- 
plications on a variety of astrophysical phenomena. In this section, we discuss 
two: the effect on the propagation of ultra-high energy cosmic rays (UHE- 
CRs) and the inducement of Faraday rotation. 



6.1 Propagation of ultra- high energy cosmic rays 

UHECRs are known to originate from extragalactic sources. Hence, on their 
path through the intergalactic space, the trajectories of UHECRs are de- 
fiected by the magnetic fields bet ween sources and us, the IGMF as well as 
the galactic magnetic field (GMFV fPas et aL I (120081^ studied the effect of the 



IGMF described in Section [5] on the propagation of UHECRs. Under the 
premise that the sources of UHECRs are strongly associated with the LSS 
of the universe, super-GZK protons of £^ > 10^^ eV were injected by AGN- 
like sources located inside clusters of galaxies. Then, the trajectories of the 
protons were followed, while taking account of the energy loss du e to inter- 
action s with the cosmic microwave background (CMB) radiation. iDas et al.l 
( 2008f ) found that the defiection of UHECR trajectories is caused mostly by 



the IGMF in filaments, rather than the the intracluster magnetic field in 
clusters; it is because filaments fill a larger fraction of volume, although their 
magnetic fields are weaker, than clusters. With the gyroradius of protons 

that is, with the gyroradius corresponding typical filament size for the mag- 
netic fields typical in filaments, the defiection due to the IGMF in filaments 
is expected to be significant. Indeed, the defiection angle between the arrival 
direction of super-GZK protons and the sky position of their actual sources 
was found to be quite large with the mean value of (9) ^ 15° and the median 
value of 7 - 10°. 

The above defiection is much larger than the deflection by the GMF; 
the deflection a ngle due to the GMF was predicted to be a few degree (see, 
e.g.. lTakami~nd Sata,200i . As a matter of fact, the deflection angle of 
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{6) ^ 15° is also much than the angular wmdow of 3.1° used by the Auger 
collaboration in the study of the c orrelation between their highest energ y 
UHECR events and nearby AGNs (iThe Pierre Auger Collaboration! 120071 ). 
Althou gh the defle ction angle is large. iRvu et al.l (|2010[ ) noticed that in the 
work of 'Das e t al.l ( 20081 ). the separation angle between the arrival direction 
of super-GZK protons and the sky position of nearest AGNs is substantially 
smaller with (S) ^ 3.5 — 4°, which is similar to the mean angular distance 
in the sky to nearest neighbors among AGNs. This mean separation angle is 
comparable to the angle used in the correlation study by the Auger collab- 
oration. This is a direct consequence of the fact that the sources and us, as 
well as the IGMF, all trace the matter distribution of the universe. That is, 
although the IGMF described in Section [S] predicts larger deflection of UHE- 
CRs, it is not necessarily inconsistent with the intervening magnetic fields 
implied in the Auger experiment. 



6.2 Faraday rotation induced by the IGMF 

The IGMF described in Section [5] induces RM. With characteristic lengths 
smaller than the dimension of clusters or filaments, the inducement of RM is 
a random walk process; the standard deviation of RM is 



O-RM = 0.81 He 5||rms ]j (^—^^ ^ ^'^d ni ^, (33) 

where n^, i^H rT psi an d L are in units of cm^"^, /iG, and pc, respectively 



(|Cho and Rvul 120091 ). i^iirms is the rms strength of line-of-sight magnetic 



field and L is the path length. Here, the coherence length for RM is given 
as I = (3/4)Lint. For clusters, with he ^ 10'^ cm^'^, Brms ~ a few /xG, 
L ^ 1 Mpc, and Lint ^ 20 kpc (see Section [51), we get q-r.m ^ 100 ra d 
m~^, which agrees with the observed RM in clusters ( Clarke et al.ll200ll) . 



The magnetic field strength in filaments quoted in Section [5] is (B) ~ 10 
nC. But the value depends on how it is averag ed; that is, {B^y/^ - a few 
X 10 nC, {pB)/{p) ~ 0.1 pG, and {{pB)Y^^ / {pY^^ - a few x 0.1 pG, in 
the warm-hot intergalactic medium (WHIM) with T = 10^ — 10^ K which 
mostly composes filaments. The average value of ((pi?)^)^/^/(/3^)^/^ should 
be relevant to RM. With fie ~ 10"^ cm^^^ ^^^^ ^ 0.3 pG, L ^ "o Mpc, and 
Lin t ^ 300 kpc, we get asM ~ 1 rad vnT^ through filaments. 

lAkahori and Rvul (|2010[ ) studied in details the RM due to the magnetic 



fields in filaments using the IGMF described in Section [SJ Figure [TU] shows 
a typical RM map in (28 /i~^Mpc)^ area; the spatial distribution of RM 
traces the large-scale distribution of matter, showing two clusters and a fil- 
amentary structure containing several groups. The resultant RM is domi- 
nantly contributed by the density peak along line of sight. The rms of RM 
through filaments at the present univer se was predicted to be ~ 1 rad vnT^, 
which agrees with the estimation above (|Cho and Rvul[2009i) . FigurefTTlshows 
the two-dimensional power spectrum of the RM in the local universe within 
100 h-^ Mpc; PRM(fc) ~ |RM(k)|2fc, where RM(k) is the Fourier transform 
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Fig. 10 RM map of path length of L = 
(28 /i^^Mpc)^ area at z = 0. 



W 

-i.m 
-5.od 

100 in the local universe o 



of RM(x). The power spectru m of RM peaks at k ^ 100 , which corresponds 
to ~ 1 Mpc. In addition, lAkahori and Rvul ()2010 f) predicted that the 
probability distribution function (PDF) of |RM| through filaments follows 
the log-normal distribution. We note that RM ^ 1 rad m~^ is an order of 
magnitude smaller than the valu es of |RM| toward the Hercules and Perseus- 
Pisces superclusters reported in IXu et"aLl (|2006l ). The diff erence is niostly 
due to the mass- weighted path length; the value quoted by IXu et all ()2006[) 
is about two orders of magnitude larger than ours. 

RM ^ 1 rad m^^ due to the IGMF in filaments is too small to be con- 
fidently observed with currently available facilities. In addition, the galactic 
foreground of ~ 10 rad m^^ (toward halo) poses an additional challenge for 
its observation. The next generation radio interferometers, however, are ex- 
pected to be able to observe the RM. Particularly, the SKA could measure 
RM for ^ 10^ polarized extragalactic sources across th e sky with an average 
spaci n g of ^ 60 arcsec bet ween lines of sight (papers in ICarilli and Rawlingsl 
l2004l ) (jKrause et all 120091 ) . enabling us to investigate the IGMF in the LSS 
of the universe. 



7 Structure Formation and Magnetic Fields 

The existence of magnetic fields in the early universe, although probably 
weak, can have consequences on the structure formation itself. In this section, 
we briefiy review recent developments in our understanding of magnetized 
structure formation. 
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Fig. 11 Two-dimensional power spectra of RM for 3 x 16 two-dimensional pro- 
jected maps from 16 simulations with different realizations of initial condition to 
compensate cosmic variance. The thick gray line shows the average. 



7.1 The linear regime 

Studies of magnetized structure formation go back to the late 1960s with the 
early efforts based on Newtonian gravity and the relativistic approaches be- 
ing a relatively recent addition. All treatments typically work within the ideal 



on density inhomogeneities (IRuzmaikina and RuzmaikinI 1971 


, Wasserman 


1978'.'Fcnnelly and Evans' 1980VKim et al.ll996l Battaner et al.l 


1997,lTsaeas and Barrowl 


1997, Tsagas and Maartcns 20003) • These effects generally come in the form 



of scalar, vector and (trace-free) tensor distortions. The former are those 
commonly referred to as density perturbations and represent over-densities 
or under-densities in the matter distribution. Vector inhomogeneities describe 
rotational (i.e. vortex-like) density perturbations. Finally, te nsor-type density 
inhomogeneities correspo nd to shape distortions Q Following ( Tsaeas and Barrowl 
Il997l iBarrow et"aIll2007^ ■ we define the scalar 



P 



(34) 



^ It should be made clear that trace-free tensor inhomogeneities (i.e. shape de- 
formations) and pure-tensor distortions (i.e. gravitational waves) are two different 
types of perturbations. 
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which describes hnear density perturbations and corresponds to the more 
famihar density contrast 5p/ p. Note that positive values for A indicate over- 
densities and negative ones under-densities. Also, a is the cosmological scale 
factor and D'^ = D"Da is the 3-D Laplacian operator that corresponds to 
an observer moving with 4- velocity WqO In a perturbed, weakly magnetized 
and spatially flat Friedmann-Robertson- Walker (FRW) universe, the above 
defined scalar evolves according to 

A^iwHA-{l + w)Z+'^cl{l + w)HB, (35) 

where over-dots denote proper-time derivatives (relative to the Ua-frame). 
Also, w = p/ p is the (constant) barotropic index of the matter, H = a/a is 
the background Hubble parameter (with — 3H there) and = / p{l+w) 
is the square of the Alfven speed. The variables Z = a^D^O and B = 
{a^ / B^)D'^ B'^ describe linear inhomogeneities in the smooth Hubble expan- 
sion and in the magnetic energy density respectively. Then, to first order, 

11 1 

Z = -2HZ --pA + - cl{l + w)pB B^A - - cl^^B (36) 



and 



2' 4 l + w 2 



B^——-A, (37) 

3(1 + w) 



respectively. Note that cj. — p/p is the square of the adiabatic sound speed 
and we have assumed that B^ p, given the relative weakness of the mag- 
netic fields. 

Equation ([55)1 shows that magnetic fields are generic sources of linear 
density perturbations. Indeed, even if A and Z are zero initially, A will 
take nonzero values solely due to the magnetic presence. Also, Equation ((57|) 
ensures that perturbations in the magnetic field energy density evolve in tune 
with their matter counterparts (i.e. B cc A). Finally, we should emphasize 
that only the pressure part of the Lorcntz force contributes to the linear 
relations ^ and 

The system of Eq uations (1551) - (1571) has an alytical solutions in the radi- 
ation and dust eras (jTsagas and Barrowlll997l ). Before equipartition, when 
w = 1/3 = c^, H = l/2t, p = 3/4:t'^ and cl = W'^/Ap = constant, large-scale 
magnetized density perturbations obey a power-law solution. In particular. 



^ For an observer with 4-velocity Ua (so that UaU°' — —1), the tensor hab = Qab + 
UaUb projects orthogonal to Ua and Da — hJ'Wb defines the covariant derivative 
operator of the spatial hypersurfaces {gab is the space-time metric and Va the 
associated covariant derivative). 

^ The Lorentz force splits as eabcB^curlB" = T>a/2B'^ - B^'DbBa, with the for- 
mer term corresponding to the magnetic pressure and the latter to the field's 
tension. The effects of the B-field on (scalar) density perturbations propagate 
via the divergence of the Lorentz force. To leading order, the latter is given by 
D" (Sabc-B'' curls'^) = D^BVS - B^TlbQ^Ba - KB'^/a^, where isT = 0,±1 repre- 
sents the background 3-curvature index. Given that D"_Ba = at the ideal-MHD 
limit, the magnetic tension effects are not included in this perturbative level, un- 
less th e spatial curvature of the FRW model is accounted for (|Tsagas and MaartensI 
I2000al lbh. 



30 



on super-horizon scale s and keeping only the dominant Krowing and deca ying 
modes, one arrives at ( Tsaeas and Maartensll2000al iBarrow et al ]|2003) 



(38) 



In the absence of magnetic fields, we recover the standard growing and de- 
caying modes of Z\ cx t and A cx respectively. So, the magnetic presence 
has reduced the growth rate of the density contrast by 4c^/9. 

Well inside the horizon we can no longer ignore the role of the pressure 
gradients. There, the fc-mode oscillates like a magneto-sonic wave with 



A 



cx sm 




(39) 



wher e Xk = a/k is the perturbed scale and Xh = \ IH the Hubble hori- 
zon ( Tsaeas and MaartensI l2000al iBarrow et al.ll2007[ ). Here, the magnetic 
pressure increases the effective sound speed and therefore the oscillation fre- 
quency. The former makes the Jeans length larger than in non-magnetized 
models. The latter brings the peaks of short-wavelength oscillations in the ra- 
diation de nsity closer, leaving a potentially observable signature in the CMB 
spectrum ( Adams et al.lll996r) . 



When dust dominates, w = Q — c], H — 2/3t, p — i/St^ and cl — B"^ / p cx 
^-2/3 fpj^gj^^ on superhoriz on scales, the main growing and decaying mo des 



t ' . iiien, on superhoriz on scales, the mam growing and decayiiig mo c 
of the density contrast are ( Tsaeas and Barrowlll997l Barrow et al.l 20071) 



A = Cit° 



(40) 



with ai^2 = ^[1 ± 5-\/l — (32/75)(ca A_f//Afc)Q]/6. In the absence of magnetic 
fields we recover again the standard solution with ai — 2/3 and a2 — — 1. 
Thus, as with the radiation era before, the magnetic presence slows down 
the growth rate of density perturbations. Also, since B (x A [see Equation 
(P?)) ]. the above describes the linear evolution of the magnetic energy-density 
perturbations as well. This means that cosmological magnetic fields trapped 
inside an overdense region of the post-recombination universe could grow by 
approximately two to three orders of magnitude. Note that the aforemen- 
tioned increase is different from the one occurring during the subsequent, 
nonlinear contraction of a protogalactic cloud (see Section 17.21 below) . 

The field pressure also leads to a magneti cally induced Jeans length, below 
which density perturbations cannot grow ( Subramanian and Barrow! Il998l 
ISethi and Subramanian|[2005l) . As a fraction of the Hubble radius, this pur ely 
magnetic Jeans scale is ( Tsaeas and Maartensll2000ai IBarrow et al.l[2007l ) 



A, 



CaA 



H 



(41) 



Setting i? ^ 10 ^ G, whi ch is the ma x imum homogeneous field strength 
allowed by the CMB (ZeldovichI [l970al IBarrow et aLlll997^ ■ we find that 
Xj ~ 10 kpc. Alternative, magnetic fields close to 10"'' G, like those found 
in galaxies and galaxy clusters, give Aj ~ 1 Mpc. The latter lies intriguingly 
close to the size of a cluster of galaxies. 
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Overall, the magnetic effect on density perturbations is rather negative. 
Although magnetic fields generate this type of distortions, they do not help 
them to grow. Instead, the magnetic presence either suppresses the growth 
rate of density perturbations, or increases the effective Jeans length and 
therefore the domain where these inhomogeneities cannot grow. This negative 
rol e of magnetic fields, which was als o observed in the Newtonian treatment 
of ( Ruzmaikina and Ruzmaikinlll97l[ ) , reflects the fact that only the pressure 



part of the Lorentz force has been incorporated into the equations. When the 
tension component (i.e. the elasticity of the field lines) is also accounted for, 
the overall magnetic effect can change and in some cases it could even reverse 
(|Tsagas and Maartengir2000bD . 



Magnetic fi elds also induce and affect rotational, vortex-li ke, density in- 
homogeneities ( Wassermanlll978LlTsagas and Maartensll2000a^ . To linear or- 



der, these are described by the vector Wa = —{a /2p)eabcOr^'^ P, with Sabc 
representing the 3-D Levi-Civita tensor. Then, on an spatially fiat FRW 
background, 

Wa = -4.HWa - ^ + J c^D^W, , (42) 
after matter-radiation equality ( Tsaeas and MaartensI [2000al iBarrow et al.l 



[2007) . Defining Aa = c^Xh as the Alfven horizon, we may write the associated 
solution in the form 

W(^k) ^ Cie' + C2t"\ (43) 



with ai,2 = — [5 ± ^1 — (48/9)(Aa/Afc)o]/6. On scales far exceeding the 
Alfven horizon, Aa ^ A^ and the perturbed mode decays as W cx t^"^/^. This 
rate is considerably slower than W oc the decay rate associated with 
magnetic-free dust cosmologies. Well inside A a, on the other hand, magne- 
tized vortices oscillate like Alfven waves, with (jTsaeas and Maarteng|[2000af l 



cos 



9 V A 



(44) 



Thus, the effect of magnetic fields on a given vortex mode is to reduce its 
standard depletion rate. Analogous is the magnetic effect on namely on 
the vorticity proper. Hence, magnetized cosmologies appear to rotate faster 
than their magnetic-free counterparts. In contrast to density perturbations, 
magnetic fields seem to favor the presence of vorticity. This qualitative dif- 
ference should probably be attributed to the fact that the tension part of the 
Lorentz force also contributes to Equation 

In addition to scalar and vector perturbations, magnetic fields also gen- 
erate and affect tensor-type inhomogeneities that describe shape-distortions 
in the density distribution (jTsaeas and Maartendl2000a[) . An initially spher- 
ically symmetric inhomogeneity, for example, will change shape due to the 
magnetically induced anisotropy. All these effects result from the Lorentz 
force. Even when the latter is removed from the system, however, mag- 
netic fields remain active. Due to its energy density and anisotropic nature, 
for example, magnetism affects both the local and long-range gravitational 
fields. The anisotropic magnetic pressure, in particular, leads to shear distor- 
tions and subsequently to gravitational- wave production f|CaDrini and DurreJ 
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120021 It saga sl l2002llWangll201fll) . Overall, magnetic fields are a very versatile 
source. They are also rather unique in nature, since they are the only known 
vector source of energy. An additional unique magnetic feature, which re- 
mains relatively unexplored, is the field's tension. When we add to all these 
the widespread presence of magnetic fields in the universe, it is not unrea- 
sonable to say that no realistic structure formation scenario should a priori 
exclude them. 



7.2 Aspects of the nonlinear regime 

The evolution of large-scale magnetic fields during the nonlinear stages of 
structure formation is addressed primarily by means of numerical simula- 
tions. The reason is the high complexity of the nonlinear MHD equations, 
which considerably hampers analytical studies, unless certain simplifying as- 
sumptions are imposed. 

The simplest approximation is to assume spherically symmetric compres- 
sion. Realistic collapse, however, is not isotropic. In fact, when magnetic fields 
are present, their generically anisotropic nature makes the need to go beyond 
spherical symmetry greater. Certain aspects of anisotropic contraction can b e 
analytically studied within the Zeldovich approximation (IZeldovichI 1 1 9 70bf) . 
The latter is based on a simple ansatz, which extrapolates to the nonlinear 
regime a well known linear result. The assumption is that the irrotational 
and acceleration- free linear motion of dust, also holds during the early non- 
linear stages of galaxy formation. This allows the analytical treatment of 
the nonlinear equations, leading to solutions that describe anisotropic (one 
dimensional) collapse and to the formation of the well-known Zeldovich pan- 
cakes. 

Suppose that magnetic fields are frozen into a highly conductive proto- 
galactic cloud that is falling into the (Newtonian) potential wells formed by 
the Cold Dark Matter (CDM) sector. Relative to the physical coordinate 
system {r"}, the motion of the fluid velocity is Ua = 3Hra + Va, where 
H ~ d/a is the Hubble parameter of the unperturbed FRW background and 
Va is the peculiar velocity of t he fluid (with a = 1, 2, 3). Then, the magnetic 
field induction equation reads ( Bruni et al.ll2003l) 



Ba = -2HBa - ^ + a^isB^ , (45) 

where — d°'Va and ctq,^ = d(p'Va) are the peculiar volume scalar and the 
peculiar shear tensor respectivelyO The former takes negative values (i.e. 
1? < 0), since we are dealing with a protogalactic cloud that has started to 



* The Newtonian theory is a very good approximation, since we are dealing with 
non-relativistic matter and the scales of interest are well inside the curvature radius 
of the universe. 

^ When dealing with purely baryonic collapse, the Zeldovich ansatz only holds 
during the early stages of the nonlinear regime, when the effects of the fluid pressure 
are negligible. Assuming that the contraction is driven by non-baryonic CDM, 
means that we can (in principle) extend the domain of the Zeldovich approximation 
further. 
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turn around and collapse. Also note that the first term in the right-hand side 
of Equation (HSl) reflects the background expansion, the second is due to the 
peculiar contraction and the last carries the anisotropic effects. Introducing 
the rescaled magnetic field Ba = a^Ba, the above expression recasts into 

B'^ = -^dBa.+a^^B^ , (46) 

with primes indicating differentiation with respect to the scale factor. Also 
?? — aHd and Oap — aH&ap, where 19 = d°'Va and Uap = di^pv^) (with 
Va = ax'^ and Va = aHva). In the shear eigen- frame, dap = (o'n, a'22, CT33) 
and Equation (|46| leads to 

B[ = -^m + i^iiBi, (47) 

with exactly analogous relations for the rest of the magnetic components. 
The resulting system describes the second-order evolution of magnetic fields, 
which is frozen-in with the highly conductive matter of a collapsing proto- 
galaxy, within the limits of the Zeldovich approximation. To obtain analytical 
solutions, we recall that in the absence of rotation and acceleration, the pe- 
culiar volume scalar is given by 

d = + + , (48) 

l + aXi l + aX2 l + aXs ^ ' 

At the same time, 

^11 = -^-^^, (49) 
1 + aXi 6 

while analogous expressions hold for the other two shear eigenvalues. Note 
that Ai, A2 and A3 are the e igenvalues of the i nitial tidal fi eld that determine 
the nature of the collapse ( Matarres3ll996l iBrunil [19961 ) . One-dimensional 



collapse along, say, the third eigen-direction is characterized by Ai = = A2 
and by A3 < 0. In that case, the pancake singularity is reached as a — >■ — I/A3. 
Spherically symmetric collapse, on the other hand, has Ai = A2 = A3 = A < 
0. Then, we have a point-like singularity when a — 1/A. 

Substituting, Equations (pS]) and ([^^ into the right-hand side of Equation 
(|T7)l we arrive at the solution 



Si = 



(l-haoA2)(l + aoA3) 



(l-|-aA2)(l + aA3) 



(50) 



for the first of the magnetic components. A similar calculation leads to exactly 
analogous equations for B2 and B^. The zero suffix in the above indicates 
a given time during the protogalactic collapse. The ratio ao/a in parenthe- 
ses refiects the magnetic dilution due to the background expansion and the 
brackets monitor the increase in the field's strength caused by the collapse 
of the protogalactic cloud. According to (l50l) . when dealing with pancake 
collapse along the third eigen-direction, the i?3-component decays as a~^, 
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while the other two increase arbitrarily as a —> — I/A3. Alternatively, during 
a spherically symmetric contraction, the magnetic field evolve as 

Here, all the magnetic components diverge as we approach the point sin- 
gularity (i.e. for a — >■ — 1/A). Comparing the two results, we deduce that 
the anisotropic (pancake) collapse leads to a stronger increase as long as 
A3 < A. The later is always satisfied, provided that the initial conditions 
are the same for both types of collapse, given that A3 = ^^o/(l — aof^o) and 
A = i?o/(3 — aoi^o) [see expression (|48|) above]. 

Our qualitative analysis indicates that magnetic fields trapped in an 
anisotropically contracting protogalactic cloud will increase beyond the limits 
of the idealized spherically symmetric collapse. Also, the amplified magnetic 
fields will end up essentially confined to the galactic plane. Quantitatively, 
the achieved final strength depends on the ti me the magnetic back-reaction 
has grown strong enough to halt the collapse ( Zeldovich et al.l 1983). Thus, 



the longer the anisotropic collapse persists, the stron ger the residual mag- 
netic fields. The analytical study of (jBruni et al.ll2003l) . in particular, showed 
that (realistically speaking) the anisotropy could add one or two orders of 
magnitude to the magnetic strength achieved through conventional isotropic 
compression. These results appear in very good agreement with numerical 
studies simulating shear and tidal effects on the magnetic field evolu tion in 
galaxies and galaxy clusters ( Roettieer et al.l[l999L iDolag et al ][2QQ2). 



8 Summary 

The cosmic web, a network of filaments and nodes wherein most galaxies 
reside, is a prediction of the highly successful vlCDM cosmology and appears 
to be borne out by by simulations and observations. The web is filled with 
an ionized plasma, the IGM, which is expected to be permeated by magnetic 
fields. In this chapter, we reviewed recent developments in our theoretical 
understanding of the nature and origin of IGMF with a special focus on 
fields outside of clusters. We addressed two basic questions: First, "Can the 
process of structure formation generate seed magnetic fields and amplify 
them?" and second, "What role do magnetic fields play in structure formation 
and in other astrophysical processes?" We presented several plasma physics 
mechanisms for the generation of seed magnetic fields, and showed that these 
fields could be amplified during the first star formation and later during the 
formation of the cosmic web. We saw that magnetic fields with a strength 
of (B) ^ 10 nG are expected in filaments at the present universe, while 
magnetic fields should be stronger in and around clusters/groups. We then 
discussed the effects and implications of magnetic fields on the formation 
of the first stars. We also presented a formal treatment of the evolution of 
density perturbations in the presence of magnetic fields and their effects on 
the formation of structures. 



35 



Magnetic fields in the IGM arc difficult to be observed, mainly due to the 
weak nature of them, as well as due to the diffuseness of the media. However, 
the next generation radio interferometers including the Square Kilometer Ar- 
ray (SKA), and upcoming SKA pathfinders, the Australian SKA Pathfinder 
(ASKAP) and the South African Karoo Array Telescope (MeerKAT), as well 
the Low Frequency Array (LOFAR) will enable us to investigate magnetic 
fields outside clusters with high-se nsitivity observations of sy nchrotron radi- 
ation and RM (see, e.g., papers in ICarilh and Rawlind[200l . 
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